ON THE VANISHING VISCOSITY LIMIT IN A DISK 



JAMES P. KELLIHER 

Abstract. Let u be a solution to the Navier-Stokes equations in the 
unit disk with no-slip boundary conditions and viscosity v > 0, and let u 
be a smooth solution to the Euler equations. We say that the vanishing 
viscosity limit holds on [0, T] if u converges to u in L°° ([0, T]; L 2 ). We 
show that a necessary and sufficient condition for the vanishing viscosity 
limit to hold is the vanishing with the viscosity of the time-space average 
of the energy of it in a boundary layer of width proportional to v due to 
the modes (eigenfunctions of the Stokes operator) whose frequencies in 
the radial or the tangential direction lie between L(u) and M{v). Here, 
L{y) must be of order less than 1/u and M{u) must be of order greater 
than \jv. 



1. Introduction 

In the presence of a boundary, the question of whether solutions of the 
Navier-Stokes equations with no-slip boundary conditions converge to a so- 
lution of the Euler equations as the viscosity vanishes — the so-called van- 
ishing viscosity limit — is very difficult. The convergence of most interest is 
of the velocities, uniformly over finite time and I? in space. Except in the 
very special case of radially symmetric initial vorticity in a disk, where con- 
vergence is known to hold (see Theorem 16. ip , the question of convergence 
or the lack thereof is unresolved for nonzero initial velocity in a bounded 
domain. (For a half-space with analytic initial data, the vanishing viscosity 
limit is shown to hold in |14j.) 

Tosio Kato in [6] gave necessary and sufficient conditions on the velocity 
u of the Navier-Stokes equations for the vanishing viscosity limit to hold. 
The most interesting of these is that 

cT 

||Vt*(<)|||,a (r 5 dt -» as v -> 0, 



where T cu is the boundary strip of width cu with c > fixed but arbitrary. 
Making only a small change to Kato's proof, it is possible to replace Vit 
with the vorticity u = uj{u) = d\u 2 — 02U 1 , giving Equation (|2.3p (see [7]). 
(The necessity of Equation (12.31) is immediate from Kato's condition, but 
because we do not have a boundary condition on the inner boundary of T cu 
the sufficiency of the condition requires proof.) 
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Other necessary and sufficient conditions were established by Teman and 
Wang in [15] and [17] . These are the conditions in Equation (12. 5ft and Equa- 
tion (12. 6p of Theorem 12.31 an d involve only the derivatives in the directions 
tangential to the boundary of either the tangential or normal components of 
the velocity, though for a slightly larger boundary layer. Finally, a condition 
that requires that the average energy density in the boundary layer of the 
same width as Kato's vanish with viscosity, Equation (j2.T|) . is proven in [7]. 
All these conditions (which apply to a bounded domain in dimensions 2 and 
higher) are summarized in Theorem 12.31 

We consider the issue of vanishing viscosity in the (unit) disk and look for 
weaker necessary and sufficient conditions for the limit to hold. The reason 
for working in the disk is that the simple geometry allows us to make quite 
explicit calculations using the eigenfunctions of the Stokes operator, which 
are composed of Bessel functions of the first kind. In a sense, this connects 
the energy method with the geometry. What we find is that we need only 
consider certain ranges of frequencies (or equivalently, length scales) in the 
various conditions: this is Theorem l2.41 Although Theorem l2.4l is specific to 
the disk, there is no hydrodynamical reason to expect the disk to be special 
as regards the vanishing viscosity limit, so one would expect a version of 
the theorem to apply to all sufficiently smooth bounded domains in M 2 , and 
probably in higher dimensions as well. We discuss this issue more fully in 
Remark (pTT|) . 

In [3], Cheng and Wang obtain a result regarding vanishing viscosity in 
two dimensions analogous to Equation (12.150 and Equation (12.160 . Their 
result applies to an approximating sequence to a solution of the Navier- 
Stokes equations as the viscosity vanishes, whereas our result applies to 
the necessary and sufficient condition for the vanishing viscosity limit to 
hold. While for the other conditions in Theorem 12.41 we use very different 
techniques than those in [3], our proof of the necessity and sufficiency of 
Equation (|2.15j) and Equation (|2.16j) uses the key inequality in their paper. 
Section [7] contains a brief comparison between the two results. 

In [13] , the authors consider the Stokes problem (linearized Navier-Stokes 
equations) external to a disk with time-varying Dirichlet boundary condi- 
tions, showing that the vanishing viscosity limit holds. In fact, they do 
much more than this, giving an explicit construction of the solution to the 
Stokes problem and showing that it can be decomposed into the sum of 
the solution to the linearized Euler equations, the solution to the associated 
Prantdl equations, and a small correction term. The symmetry of the ge- 
ometry allows the authors of [13] to construct the solutions in an explicit 
form (involving Bessel functions of the first and second kind) . The nonlinear 
term in the Navier-Stokes equations makes an explicit solution impossible 
for us; however, we can expand the solution in terms of eigenfunctions of 
the Stokes operator for which we have an explicit form (in terms of Bessel 
functions of the first kind) which we can use to obtain finer estimates on the 
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behavior of the Navier-Stokes equations in the boundary layer than would 
be possible for a general domain. 

A word on notation: We use C to represent an unspecified constant that 
always has the same value on both sides of an equality but may have a 
different value on each side of an inequality. 

2. Definitions and Kato-type conditions 

We now give definitions of the Euler and Navier-Stokes equations, and state 
the results from [6], [7J, [IS], and [T7] that we will need. 

In Section 2] we will specialize to the unit disk, but for now we assume 
only that Q is a bounded domain in M 2 with C 2 -boundary T, and we let n 
be the outward normal vector to T. 

A classical solution (u, p) to the Euler equations satisfies, for fixed T > 0, 

, p , ( d t u + u ■ Vu + Vp = f and div u = on [0, T] x $7, 
^ ' \ u- n = on [0,T] x T, and u = v° on {0} x 0, 

where divu = 0. These equations describe the motion of an incompressible 
fluid of constant density and zero viscosity. 

We assume that w° is in C k+€ {n), e > 0, and that / is in C k ([0, t] x Q) for 
all t > 0, where k = 1 or 2. Then as shown in [8] (Theorem 1 and the remarks 
on p. 508-509), there exists a unique solution u in C^ oc ([0, oo); C k+e (Q)). 

The Navier-Stokes equations describe the motion of an incompressible 
fluid of constant density and positive viscosity v. A classical solution to the 
Navier-Stokes equations can be defined in analogy with (E) by 

J dtu + u ■ Vm + Vp = vAu + f and div u = on [0, T] x Q, 
\ u = on [0, T] x T, and u = u° on {0} x O. 

We will work, however, with weak solutions to the Navier-Stokes equations. 

Definition 2.1 (Weak Navier-Stokes Solutions). Given T > 0, viscosity 
v > 0, and initial velocity < in H, u in L 2 ([0, T]; V) with d t u in L 2 ([0, T]; V) 
is a weak solution to the Navier-Stokes equations if -u(0) = u® and 

{NS) / d t u ■ v + / (u- Vu) -v + u I Vu-Vv= / fv 
Jn Jn Jn Jn 

for all v in V. (The spaces H and V are defined in Section O) 
Definition 2.2. We say that the vanishing viscosity limit holds if 

u -> u in L°°([0,T];L 2 (O)) as v -> 0. (2.1) 

Theorem 12 . 31 applies to a bounded domain with C 2 -boundary in d > 2. 
The conditions in Equation (|2.2j) and Equation (|2.4B are due to Kato ([6]), 
the conditions in Equation (j2.3f) and Equation (|2.7p appear in [7], and the 
conditions in Equation (12.51) and Equation (|2.6p are due to Temam and 
Wang ([E], [T7J). 
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Theorem 2.3. Let T > and assume that is in H and that vP is in 
C k+€ (Q), e > with k = 1 or 2. In addition, assume that 

(a) u% -» u° in L 2 (Q) asu^O, 

(b) / is in L X ([Q,T\\I?(Q)), 

(c) 11/ - /IIli([o,t] ; l2(q)) -» as 1/ -» 0. 

Lei 5 : [0, oo) — > [0, oo) 6e suc/i i/tai <5(f) converges to while 5{v)/v diverges 
to oo as — > 0. Then the the vanishing viscosity limit (Definition \2.2\) holds 
if and only if any of the following conditions holds: 

rT 

\Hs)\\ 2 LHn) ds^0 asu^O, (2.2) 
ll w ( s )lli 2 (r c „) ds -> as v -> 0. ( 2 - 3 ) 

T 

l|V^)||| 2(rci/) ds^O (2.4) 
1/ / ||V T -u T (s)||^ 2(ri(j/)) ds -> as -» 0, (2.5) 

•/ 

^ / ||VrUn(s)||i2 (r(5(iy) ) ds -> as i/ -> 0. (2.6) 

idere V7 represents the derivatives in the boundary layer in the directions 
tangential to the boundary, u T is the projection of u in the direction tangen- 
tial to the boundary, and u n is the projection of u in the direction normal 
to the boundary. 

When k = 2, these conditions are also equivalent to 

T 

W u ( s )\\l*(t,,a ds -> as v -> °- ( 2 - 7 ) 



The quantity in Equation (|2.T|) is proportional to the space-time average 
of the energy in the boundary layer. 

We show (see Remark (]5.2p ) that in Equation (12. 2p . Equation (|2.4|) . and 
Equation (12. 7\i . contributions from the high frequency modes can be ignored. 
This result applies to an arbitrary bounded domain in M. d , d > 2, with a 
C 2 -boundary. 

Our main result is Theorem 12.41 which is an improvement of Theorem 12.31 
in the special case of the unit disk. In what follows we decompose the 
solution u in the form 

00 00 

u{t,x) = ^2g m j(t)u mj (x), 

m=0 j=l 
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where (u m j) are the eigenfunctions of the Stokes operator described in Sec- 
tion [3] and Section [H and let 

N N 

u N {t,x) = Y,Smj(t)u mj (x) (2.8) 

m=0 j=l 



and 

u N 



N oo 

(M) = Y ^2gmj{t)u mj (x) (2.9) 

771=0 j = l 

with vorticities u, > (t, x) = uj(u N (t, x)) and u (t, x) = u(u (t, x)). 

As we will see in Section HI the frequency of u m k in the tangential direction 
is m and the radial frequency of u m k is, in effect, k. Thus, u N includes the 
contributions from all modes with both frequencies less than N, while u N 
includes the contributions from all modes with tangential frequency less than 
N. 

Theorem 2.4. Assume that O is the unit disk and make the same assump- 
tions on the initial data, forcing, and the function 5 as in Theorem \2.SX Let 
L and M be any functions mapping (0, oo) to Z + with 

vL(v) -> 0, vM(v) -> oo as v -» 0. (2.10) 

Then the the vanishing viscosity limit (De finition \2.2\) holds if and only 
if any of the following conditions holds: 

v f T |Ks) M(l/) - uj L{u \s)\\ 2 L2{n) ds -> as v -> 0, (2.11) 



T 

\\lo{s) - u) L{u) {s)\\ 2 L2{n) ds -> as v -> 0, 
v J \\u(s) - u L M(s)\\ 2 L2{rcu) ds^Oasu^O, 
v \\Vu M{ - v \s) - Vu Li ~ v) {s)\\ 2 L 2 {Tci/) ds ^ asu^O, 
v I ||V t «t(s) - ^tu^ & \s)\\ 2 l2( y 5(v)) ds -> as v ->• 0, 

J 
rT 

v \ \\V T u n (s) - V T 5n \ s )\\L 2 (r s(l/) ) ds —>0 as v — > 0. 
When k = 2, these conditions are also equivalent to 

- ! T \\u M ^{s) - u L ^(s)f LHTci/) ds^Oasu^O. 



2.12) 
2.13) 
2.14) 
2.15) 
2.16) 

2.17) 



6 



JAMES P. KELLIHER 



Observe, for instance, that u M ^ —u L ^ in Equation (|2.17p represents the 
contribution from all modes whose frequencies in the radial or the tangential 
direction lie between L{y) and M[y). 

Remark 2.1. By Lemma IA.3I and Equation (|4.7p . u N is essentially the 
contributions of all the modes with eigenvalues less than CN 2 . In fact, 
suppose that we replace the definition of u N in Equation (|2,8p with 

u N (t,x)= Yl 9j(t) Uj (x), (2.18) 

the single subscripts in Equation (|2.18p referring to the eigenfunctions and 
eigenvalues of the Stokes operator on a general domain in M, d , d > 2, defined 
in Section [3J It follows easily from Theorem 12.41 that the conditions in Equa- 
tion (|2.1ip . Equation (|2.13p . Equation (|2.14|) . and Equation (|2.17p continue 
to be equivalent to the vanishing viscosity limit. It is in this form that we 
would expect Theorem 12.41 to generalize to fairly arbitrary smooth domains 
in M 2 and — with N 2 in Equation (|2.18p replaced by N raised to some other 
power — to domains in R rf , d > 3. The obstacle to establishing this gen- 
eralization is the difficulty of obtaining the equivalents of Lemma IA.8I and 
Lemma I A. 91 — along with an approximate form of Lemma IA. 101 — for high 
frequencies. 

3. The Stokes operator in a bounded domain 

Before specializing to the case of a disk, we discuss first some general prop- 
erties related to the Stokes operator. 

We define the function spaces H and V as follows (see Section 1.1.4 of 
[16j for more details). First let 

V={ue (D(ft)) 2 : divu = 0} 

be the space of vector- valued divergence- free distributions on ft. We let H 
be the closure of V in L 2 (ft) and V be the closure of V in Hq ($7). Alternate 
characterizations of H and V are 

H = [u G (L 2 (Q)) 2 : divu = in ft, u • n = on T} , 

V = {u e (L 2 (n)) 2 : divu = in ft, u = on V) , 

the boundary conditions applying in terms of a trace. 

By (•, •) we mean the inner product in L 2 (ft): (/, g) = f n fg. (It will be 
convenient to use complex-valued eigenfunctions, so the complex conjugate 
is required in this definition. Our velocity fields and vorticities, however, 
are real, so conjugation will not always appear in our calculations.) Then 
(u,v) H = (u,v) and (u,v) v = (X7u,Vv). 

Although V is dense in H it is not dense in H n H 1 (ft) (with the H 1 - 
norm). Informally, this is because each element of V is zero on T and so the 
limit of a sequence of elements in V cannot become nonzero on the boundary 
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without the gradient near the boundary becoming indefinitely large. More 
formally, we have Lemma I3.1L 



Lemma 3.1. The space V is not dense in H n H 1 ^!). 

Proof. Let u be any element of V. Then its vorticity u is in L 2 (J7) C L l (Q) 
and must satisfy 



oj= I Ail>= / div = / VV> • n = - / u • n = 0, (3.1) 
n Jn Jn Jr Jr 

where ip is the stream function: u = V^V = ( — ^2'0)^1'0) an d w = 
Because only u-rv = for it in .ff, the same cannot be said for it in HnH 1 ^): 
Let it in HnH 1 ^) have vorticity a; with nonzero total mass. Then for any 
sequence {vj} in V, 



l^-^i)ll L 2 (n) >C||o;-a;(^)llii ( n) > c 



w — U)(Vj)) 




C 



>0, 



so V cannot be dense in H n -fT 1 (ri) . □ 

We now briefly describe the properties we will need of the Stokes operator 
A on f2, referring the reader, for instance, to Section 1.2 of [16J for more 
details. One way to define A is that given u in V P\H 2 (Q), Au in H satisfies 
Au = —Au + Vp for some harmonic scalar field p. We have D(A) = V n 
H 2 (Q) with vl mapping D(A) onto -ff , and there exists a set of eigenfunctions 
{uj} for j4, complete in H and in V, with corresponding eigenvalues {A^}, 
< Ai < A2 < • • • , and each Uj is in H 2 (£l) since we are assuming that T is 
C 2 . (When we specialize to the disk, the eigenfunctions will be in C°°(0).) 
An eigenfunction Uj of A satisfies Auj = XjUj or, equivalently, 

Auj + XjUj = Vpj, Apj = 0, div Uj = 0, on f£, - , 

Uj = on r. \ ■ J 

The eigenfunctions are orthogonal in both H and V . The usual convention 
is to make the eigenvectors orthonormal in H, but we will find it more 
convenient to normalize them to be orthonormal in V so that IIVttjH^^ = 

HWl^n) = 1 and 

2 111 

Imh^u) = ( u ji u i) = 7- ( u J' Au j) = 7- ( Vn i' Vn i) = T~- ( 3 - 3 ) 

Moreover, we have Lemma 



Lemma 3.2. If u is in V with u = uo{u) then 

00 

u = J2{ UJ ,Uj)u j , (3.4) 
3=1 

with the sum converging in both V and H . 
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Proof. Let u be in V and let u n = Y?j=i({ u ) u j) h I ( u j' u j) n) u y Then u n 
converges in H to u because {%} is complete in H. But, 

(u,Uj) H Xj(u,Uj) {u,Auj) (Vu,S7uj) 
(uj ,Uj) H Aj {uj , uj ) (uj, Auj ) ( Vuj , Vuj ) 
= (Vu, Vuy) = (uj,ujj) , 

so the expansion of u in V in terms of the eigenfunctions of A is the same as 
the expansion of u in H (and the coefficients are as given in Equation (13, 4j) ). 
meaning that u n converges in V to u as well. □ 

In the proof of Lemma 13.21 we used the identity (Vu, Vt;} = (uj(u),uj(v)) 
for all u, v in V, which follows by integrating by parts. Were we to use 
the definition of u(u) as the antisymmetric matrix (Vu — (Vu) T )/2, which 
is usual in higher dimensions, this would have introduced a factor of 2 into 
Equation ([31| . 



Corollary 3.3. If u is in V then 

oo oo 

Vu = {to, ujj) Vuj and oo = (uj, luj) ojj, 
3=1 3=1 

with the sums converging in L 2 (£2). 

Since the solution u to (NS) lies in V for all positive time, we can write 



w(t) = ^gjitpj, u{t) = Yj9j(.t)i 

3=1 3=1 



A, 



(3.5) 



j=l j=l ■ 



where g 3 - are functions of time. The expansion of u will converge for all t > 
and that of u; for t > — and also for t = if and only if the initial velocity 
is in V; in general, we only assume that it is H. Because u(t) — > u° in L 2 (£l) 
as t — > 0, each gj(t) is continuous at t = 0, though this does not mean that 
u(t) is continuous in L 2 (Q) at t = 0. Also, note that 5j(t) is complex- valued 
since the eigenvectors are complex-valued, but u(t) and u(t) are real-valued. 

4. Eigenfunctions of the Stokes operator in the unit disk 

We now fix to be the unit disk in M. 2 centered at the origin. 

In [TU] , a complete set of eigenfunctions for the annulus is derived in terms 
of Bessel functions of the first and second kind, J n and Y n . By ignoring the 
terms involving Y n and modifying somewhat the calculation of the eigenval- 
ues, one can easily obtain the eigenfunctions for a disk. We will, however, 
derive the vorticity of the eigenfunctions directly, as this is quite easy. In 
determining the eigenvalues and the velocity of the eigenfunctions, which is 
more difficult, we will rely on the results in |10j . 
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Taking the curl of Equation (13. 2D (with u = Uj), we see that the vorticity 
uj = u>(u) satisfies 

Au + Aw = on O, , . 

u = on r. \ ■ J 

That is, uj is an eigenfunction of the negative Laplacian, but with boundary 
conditions on the velocity u. 

Ignoring for the moment the issue of boundary conditions, we use sepa- 
ration of variables to look for a complete set of solutions to Au; + Xuj = 
on £1. Writing 

u(r,9) = Ur)e ind 

in polar coordinates, n = 0,1,2, ... , and using 

d 2 Id Id 2 
A = — r + -— + 



dr 2 r dr r 2 d6 2 
gives 



d2fn +-^+(x-^)f n )e^ = 0. (4.2) 



dr 2 r dr \ r 2 
Since J n , the Bessel function of the first kind of order n, is a solution of 

d 2 Us) + lJJM( l _^Us) = 0, (4.3) 



ds 2 s ds \ s 2 

making the change of variables s = A 1//2 r, we see that Equation (|4.2|) holds 
with f n (r) = J n (A x / 2 r). Thus, the eigenfunctions have vorticity of the form 
J n (\ l / 2 r)e ind and it remains to determine the eigenvalues that satisfy u = 
on the boundary. 

The easiest way to do this is to use the expressions in |10j . For n = 0, we 
drop the term involving Y\ in Equation (30) p. 406 of [10] . giving 

u k(r,9) = A ofe /2j i( A ofc 2r )^' 

1/2 

where A Q ^, , k = 1,2, ... , are the eigenvalues described below. For n > 1, 
dropping the terms involving the Bessel functions of the second kind from 
the last equation on p. (406) of [ID] , we have 

u nk (r,0) = { ^UX^r) + ^! r -A f% 

A nk J 

2 \ 

n— 1 I „in0' 




J„ +l (A^ r )-J„_ l( A^))-^»-' e 



A nk j 

the eigenvalues \ n k, k = 1, 2, ... , being described below and the D n k being 
undetermined constants. In both cases we scaled the eigenfunctions differ- 
ently than in [10J. A direct calculation shows that 
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where C nk is a normalization constant. A direct calculation also shows that 
div u nk = 0. 

For n = 0, A^ 2 = j lk , where 

j n k is the k-th positive root of J n+ \{x) = 0, (4-4) 

as this gives UQ k (l, 0) = 0. Setting u nk (l, 0) = we obtain two equations in 
the two unknowns X nk and D nk . We eliminate D nk from the two equations to 
obtain a single equation for X nk . Then using the identity in Equation (|A.2|) 
we obtain the equation 

^ 2 4(Ay fc 2 )-nJ„(A^ 2 ) = 0. 

1/2 

Thus, A fe is the fc-th positive root of 

xJ' n (x) - nJ n (x) = -x,J n+ i(x) = 0, (4.5) 

/c = 1,2,..., where we used Equation (|A.4|) . That is, A fe = jn+l,fc- It 
follows then that 



Anfc>^n(A n ^ 



n n 
Since we are normalizing the eigenfunctions so that (uj m j,uj nk ) = 5 mn 5j k , 
we must choose C nk so that 

«i 

Jn6\\2 

n+l,k> )C 

r 2 



C nk = \\JnUn+i,kr)e my \\ 2 L2{n) = 2ir J rJ n (j n+1;k ry dr 

1 

= 2lT— [J n (jn+l,kr) 2 ~ Jn-l{jn+l,kr)Jn+l{jn+l,kr)] 



= ^Jn{jn+l,k) 2 - 

Here we used Equation (|A,9|) , 

To summarize, the vorticity of the eigenfunctions is given by 

LO nk (r,9) = C nk Jn(jn+i,kr)e md , 

with eigenvalue 

Anfc = in+l,fci (4-7) 

and where 

C nk = 1 

n = 0,1,..., k = 1,2,.... With our choice of normalization of the eigen- 
functions (Equation (|3.3p ). the velocity becomes 

, m Jn{ar) - J n {a)r n . in6 ^ 
Unk(r,0) = — , — me e r 

ir i/2 a 2 |J„(a)|r 

a ( J n+ i{ar) - J ra _i(ar)) + 2nJ w (a)r n ~ 1 ine _ 
+ 2vri/2«2| Jri ( a )| e e *' 



where a = j n+ i ik - 
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5. Proof of Theorem 12.41 
From the fundamental energy equality for (NS) we have for all t in [0, T] , 

v J l|V«||ia (n) = v j \\v\\h(n) ^ ^\\ u I\\h + 4 ll/llii([o,T];L2(n)) • 

It follows from Equation f)3.5j) and assumptions (a) and (b) of Theorem 12.31 
that for all sufficiently small v > 0, 

v I H^llism) = ^ / J2 J2 \9mj(s)\ 2 ds < C. (5.1) 

J ° J ° m=0j=l 

Theorem 5.1. With the assumptions of Theorem \2.4\ 

ton* f 11^11^0^=0. (5.2) 



Proof. Using Lemma lA.10| 



" I H^lli-cr. 



) 

t L(u) L(u) L(u) L(u) 

= V j YYYY 9mj(s)gnk(s) ds (u m j 
J ® m=0 j=l n=0 k=l 

,L{u)L{u)L(u) 

= V j YYY 9nj(s)gnk(s) ds (uj nj ,uj nk ) L2{Tcv) 
Jo n=0 j=l k=l 
t L{v) L(u) L(v) 

< v j Y Y Y \9nj(s)\ \g nk (s)\ ds \\uJ n j\\ L 2 ircv) ||wn*||xa(r. 

Jq n=0 j=l k=l 

,L(u) (L{u) \ 2 

= V Y \Y \9nj(s)\ \\Unj\\ L 2 {Tcv) ds 
J ° n=0 \j=l j 

t L(u)L(v) L(v) 

<»l YYMs)fY 



'° n=Q j=l j=l 

where we used the Cauchy-Schwarz inequality in the last step. 
By Lemma lA.3|. 

l/L{u) < C/(L(u) + 2) < C/j L{u)+1;1 = CA-Jg r 

Since vL(y) — > as v — > 0, for all sufficiently small v we have cv < C ' /L(y) < 

27rX L(!)i ^ 27rX j,i /2 for a11 j - So b y Lemma EH 

L(u) 

Y IMli 2 (r c „) ^ CvL{y). (5.3) 
i=i 
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Then using Equation (15. ip . 

L{u)L(u) 



7 

Jo 



^ W ll| 2( r c „) < CWW J EE t3n,(s)| 2 ^ < CW^), 

\ JO n — n o— 1 / 



n=0 j=l 

which vanishes with v by the condition in Equation (I2.10p . and Equa- 
tion (|5.2|) therefore holds. □ 

Remark 5.1. We could try to improve Theorem 15.11 by using u{u ) of 
Equation (|2.9p in place of ui N , thereby incorporating all of the frequencies 
in the radial direction for a given angular frequency. Unfortunately, the 

1 1 2 

best bound that one can achieve on H^VjIl^r,,) f° r j > n is the extension 
of Lemma lA.8l described in Remark (jA.ip . and this is very much insufficient 
to bound the terms with j > n. 

Another possible approach is to try to incorporate the destructive interfer- 
ence that occurs in the inner product of two eigenfunctions in the boundary 
layer that the use of Holder's inequality in our proof of Theorem 15.11 ignored. 
The best bound one can hope to obtain is that 

C 

I {^nj,^nk) L ^(r s ) I - \k-j\ 

for all 5 in [0, 1] and without restriction on n, j, or k except that j ^ k. 
We could then follow the obvious approach of decomposing the equivalent of 
the first sum in the proof of Theorem 15.11 into four pieces: a diagonal term 
where m = n and j = k and three terms containing low frequencies in j and 
k, low frequencies in j and high frequencies in k, and high frequencies in 
both j and k. If we do this, however, we will find that the factor of 1/ \ k — j\ 
is just insufficient to obtain convergence. 

Corollary 5.2. The conditions in Equation \2. 1\) and Equation h2.13\) of 
Theorem \2.J\ are equivalent. 

Proof. That Equation (|2.ip implies Equation (|2,13p follows directly from 
Theorem 12.31 So assume that Equation (|2.13p holds. Because + < 
2 ||^4|| +2||B|| for any norm, 

uj* IMll^ < 2* jf* \\u L M\\l HVm) + 2vJ* \\v-v L M\\h<r„y 

This vanishes with v by Theorem 15.11 and Equation (|2.13p , showing that 
Equation (|2.3p holds and hence by Theorem 12.31 that Equation (|2.ip holds. 

□ 



Theorem 5.3. With the assumptions of Theorem \2.4 



limi/ \\u( S )-u M ^( S )\\ 2 LHr ^ds = (5.4) 

J 
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and 

Km- f\\u Liu \s)\\* ) ds = 0. (5.5) 



Proof. We can write u(t) - u M{ - u \t) = A(t) + B(t), where 

M(u) oo oo oo 

A(t) = ^2 9mj(t)u m j(x), B(t) = ^2 ^2g m j(t)u mj (a 

m=l j=M (u)+l m=M (v)+l j=l 

and 

Ht) _ u M M{t)\\W„) < 2 \\A(t)f LHVcu) + 2 \\B(t)f LHVcu) . 

Now, 

M (f) oo 



\\A(t)\\ 2 LHTci/) < \\A(t)\\ 2 L2{u) =J2 E H" 



2 n || 2 

m illL 2 (n) 

m=l j=M(f)+l 



M(v) oo | , M 2 1 M(y) oo 

E E ^f^r^E E 



< j-^— |K*) -w MM (0[li» ( 

where we used Equation (|3,3p . Similarly. 

A 



By Equation (|4.7j) and Lemma lA.31 A^(^)i and Ai^(^) are both bounded 
below (and above) by CM(v) 2 , so 



Then, 



Ht)-u M ^{t)\\l, {Vov) < ^\\"{t)-" M{v \mh { ny 

l[\\u{s)-u M ^{s)\\l Hrcv) ds 

\u(s)-co M M(s)\\l 2m ds 



c 

< 



< 



vM{v) 2 

( ;v I \\u(s)\\h(n) ds ^ 



where in the last inequality we used Equation (|5.ip . This vanishes with v 
by the assumption on M in Equation (|2.10|) giving Equation (|5.4|) . 
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Arguing as in the proof of Theorem 15.11 



*"illL 2 (r CI/ ) ^ s 



n=0 j=l 



< 



^ ■' n=0 j=l 



for all sufficiently small v. In the second inequality we used Lemma [A.9l and 
in the last inequality we used Equation (|5.1|) . This integral also vanishes 
with v by the assumption on L in Equation (I2.10p giving Equation (|5.5I) . □ 

Corollary 5.4. The conditions in Equation \2.1\) , Equation 112. Equa- 
tion (2. 14\ ), and Equation ( 2.17\ ) of Theorem \2.4\ are equivalent. 



Proof. For sufficiently large v, L(v) < M(v), and we have 

+ 3||n L M( S )||i 2(rcy) + 3||«(*) - u M M{s)f L * {Vcv) . 
It follows from Theorem 15.31 that 

limsup- f \\u(a)f Ti(r ^<31imsup- f \\u M{v) - u L ^\\ 2 T2(v v 

In particular, the first limsup is zero if and only if the second limsup is 
zero (the reverse inequality without the factor of 3 being trivial). Then 
Equation (|2.7p of Theorem [2T3] shows that Equation (|2,17p holds if and only if 
Equation (|2.1j) holds. The sufficiency of Equation (12. lip and Equation (|2.14|> 
for Equation (|2.ip to hold then follows from Poincare's inequality in the form 

h M(u ){s) _ u m {s) f L2{Tci/) < ^ ||Vn M H(s) _ Vu w {s) f La{Fcv) 

<Cu 2 \\Vu M ^(s)-Vu L ^(s)\\l 2{n) 

= cA\^ M{u \s)-^{s)\\l 2m . 

The necessity of Equation (|2.1ip and Equation (|2.14p follow immediately 
from Theorem 12.31 □ 



Remark 5.2. If we replace the definition of u N in Equation (|2.8p with 
that in Equation (|2.18p . then it is clear that Equation (|5.4|) continues to 
hold in any bounded domain in R 2 with a C 2 -boundary. It follows as in 
Corollary 15.41 that the vanishing viscosity limit of Definition 12.21 holds if 
and only if the condition in Equation (|2.1ip . Equation (|2.14p . or (when 
k = 2) Equation (|2.17|) holds with the term involving each of these 

conditions removed. A similar result would hold in any dimension for an 
arbitrary bounded domain with a C 2 -boundary. 
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Theorem 5.5. With the assumptions of Theorem \2.4\ 

r-T 



and 



f 

v / \\V T UT {S \s)\\h { r SM ) ds -»■ as v -> ( 5 - 6 ) 

J 

r-T 

v J ||V T 5n (<5) (s)||i2 (r4(iy)) -» as i/ -> 0. (5.7) 

Proof. In the unit disk, « T = u 61 and Vr = 5 CT , where a is arc length along 
the circle of radius r, in which r is held constant. Thus, 

„ du e 1 du e 

V T u T = -TT- = — 

and for any positive integer N it follows from Poincare's inequality that 



|VT5?(»)ll£»(r,) 



< 



1 



< 



C<5 2 



5 2 



drdO 



LHFs) (! - 5 ) 



< 



(1-^ 



t; 2 



2 

£ 2 (r*) 

2 

L 2 (H) 



But, 



q2 N oo r,2 

m=0 j=l 



drdO 



N oo ^ 

= i ^ m^2g m j(s)—u m j(r,6), 

m=0 j=l 

the last equality following from the simple dependence of u m j on 6 in Equa- 
tion gS}. Thus, 

2 



Z9 2 

,~ N 



(^(s)) 6 



< 



L 2 (Q) 



N oo 



^ m'y]g m j(s)Vu m j(r,9) 



L 2 (Q) 



m=0 j=l 

N oo N oo 

= E ™ 2 E iw*)i 2 <^ 2 EE iw*)i 2 ^ ^iwiiw 

m=0 j'=l m=0 j=l 

where we used the orthonormality of the eigenfunctions in V. 
Combining these two inequalities gives 

CN 2 S 2 



|Vx^(s)" 2 



T' U T \ s ) \\L 2 (Tg) < (l_S)2 

Then using Equation (|5.1|) . 



Vd 



f T nx7 „l(S) ( , ||2 . C7Lffl 2 ^ 2 f T 



\L 2 (n) 



CL(5) 2 5 2 



< 
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This vanishes with 8 by the assumption on L in Equation (I2.10p and hence 
vanishes with v since b vanishes with v, giving Equation (15. 6j) . The proof of 
Equation (|5.7j) is entirely analogous. □ 

The technique used in the proof of Theorem 15.51 comes from the key 
inequality following Equation (3.21) in [3]. 

Corollary 5.6. The conditions in Equation \2. Equation Ii2.12\) . Equa- 
tion t2.15\) . and Equation A2.16\) of Theorem \2.J\ are equivalent. 

Proof. This corollary can be proved much along the lines of the proofs of 
Corollary 1 5 . 2 1 and Corollary 15.41 (It is here that we use the assumption that 
8{y)/v diverges to oo as u — > 0, which is needed in applying Theorem 12.31 ) 

□ 

Together, Corollary 15.21 Corollary 15.41 and Corollary 15.61 establish Theo- 
rem [274J 

6. Radially symmetric initial vorticity 

It follows from Lemma [3. II that if an initial velocity, no matter how smooth, 
lies in H but not in V and has a vorticity whose total mass is nonzero, then 
the velocity of the corresponding solution to (NS) will be discontinuous in 
H 1 ^) at time zero. In the same way, the vanishing viscosity limit of the 
vorticity cannot hold in L°°([0, T]; L 2 (Q)) for such an initial velocity, since 
the total mass of the vorticity for the solution to (E) is conserved over time. 
This means that we might expect a different character to the vanishing 
viscosity limit when the initial vorticity has zero total mass versus when it 
has nonzero total mass. 

Indeed, this is what happens in the special case of radially symmetric ini- 
tial vorticity where, when the initial velocity is in V (which is equivalent for 
radially symmetric vorticity to the total mass of the vorticity being zero) we 
obtain convergence in L°°([0, T]; L 2 (Q)) of both the velocity and the vortic- 
ity (see [1]), whereas for initial velocity in H we obtain only the convergence 
of the velocity in this space, as in Theorem 16.11 Such convergence follows 
immediately from the conditions in Equation (|2.12p . Equation (|2.15p . or 
Equation (|2. 16j) of Theorem 12.41 The convergence also follows from the suf- 
ficiency of the conditions in Equation (|2.5p and Equation (|2.6[) as established 
in [15], since both conditions are satisfied (the gradients in the tangential di- 
rection being zero) as pointed out in [T7]. When the forcing is zero, however, 
there is a simple proof that uses only Kato's original conditions. 

Theorem 6.1. Assume that u® and u° are as in Theorem \2.3\ with (for 
simplicity) = TiP , that f = f = 0, and that u = lu(u°) is radially 
symmetric. Then the vanishing viscosity limit of Equation \2. 1\) holds. 

Proof. Because u° is radially symmetric, oj remains radially symmetric for all 
time, so lo(u-Vu) = u-Vlo = 0. Then because Q is simply connected, u-Vu = 
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for some scalar field q, and the nonlinear term in (NS) disappears. Thus, 
(NS) reduces to u u (0) = ifi and 



d t u v ■ v + v \ Vu u -Vv = (6.1) 
n Jn 

for all v in V. This is the heat equation in weak form, which is invariant 
under the transformation (v, t,x) i— > (l,ut,x). That is, if u\ is a solution 
to Equation (|6.ip with v = 1, then u u (t,x) = u\(vt,x) is a solution to 
Equation (|6.ip because 1^(0) = -ui(O) = u° and 

— ui(vt, x) ■ v (x) dx + v l VuAvt, x) ■ Vv(x) dx 
not Jn 

= v \ (dtu\){yt, x) ■ v(x) dx + / Vui{vt, x) ■ Vv{x) dx =0. 
Jn Jn 

It follows that 

/"* 2 /"* 2 r* 2 

^ y o IMs)||#(n) ds = v \\ u ^ vs )\\v{n) ds = J o W^ii^h^n) dr - 

This vanishes as — > by the continuity of the integral, because tii is in 
L 2 ([0, T]; V). The limit in Equation (|2,ip then follows from the condition in 
Equation ([221) of Theorem □ 

The proof of of Theorem 16. 11 does not yield a bound on the rate of conver- 
gence in Equation (|2.1|) . Also, without assuming that the initial vorticity is 
radially symmetric, the argument in the proof of Theorem l6.1l can be applied 
to solutions to the Stokes problem (the linearized Navier-Stokes equations) 
to show that they converge in the vanishing viscosity limit to a solution 
to the linearized Euler equations (which is just the steady state solution 
u = u°). This would be more interesting, though, if time- varying Dirichlet 
boundary conditions, for instance, could be incorporated, as in [13]. 

In the simpler case of also lying in V, the solution to (E), which 
is steady state, is zero on the boundary. This eliminates the trouble- 
some boundary term that appears in the direct energy argument bounding 
\\u(t) — u(t)\\ L 2(Q-\, giving an extremely simple proof of Equation (|2.ip . We 
give, however, a longer proof of convergence using Theorem 12.31 because it 
suggests how we might treat more general initial velocities. For convenience, 
we assume zero forcing and more regularity on the initial velocity than is 
strictly necessary. 

Theorem 6.2. Assume thatu® is radially symmetric, u® =vP is in H 3 (Q)n 
V , and there is no forcing. Then the vanishing viscosity limit of Equa- 
tion \2. 1\) holds. 



Proof. Because v° is in V n H 3 (Q), u is in ££.([0, oo); V n H 3 (n)) (see, for 
instance, Theorem III. 3. 6, Remark III. 3. 7, and Theorem III. 3. 10 of |16j). 

As observed in the proof of Theorem 16.11 u ■ Vu = Vg for some scalar 
field q. Then dtu + V(p + q) = vAu and taking the divergence of both sides 
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we conclude that A(p + q) = 0. Because of the radial symmetry, however, 
p + q is constant on T and hence is constant on O. Thus, V(p + = 
and dtu = vAu. But u = on T so <9fU = on T, and it follows that 
Ait = v~ l dtu is in V. 

Because u is in V at time zero, we can use the expansion for uj in Equa- 
tion (13. 5p for all time, including for time zero. Because uj remains radially 
symmetric over time, the expansion reduces to 

oo oo 

uj(t,r,6) = ^ 90k(t)u k(r) = ^2g ok (t)C ok J (j lk r). (6.2) 

k=l k=l 

Since Au is in V, it follows from Corollary 13.31 that Au has an expansion 
like that of Equation (|6.2p : 



Au(t, r,6) = Y J hok{t)C kMhkr). (6.3) 

k=l 

(The analogous expansion of Alo including all the eigenvectors {ui n k} fails 
to converge at t = for non-radially symmetric solutions because Am is not, 
in general, in V .) 

Since AJ (j lk r) = -XokMjikr), it follows that h ok (t) = -X k{t)gok(t) 
and then by Equation (|6.ip in strong vorticity form, 

d t uj = vAuj, (6.4) 
that g' Qk (t) = -iy\ ok g ok (t). Thus, g ok (t) = gok(Q)e- uXokt and 

oo 

(t,r,9) =J2 C °k9ok(0)e- /Xokt Jo(jikr). (6.5) 



fc=i 



But then 



\UJ 



= Ec'oWlo)^- 2 ^*!!^^^)!!!,^ 



fc=i 



<E5ofc(0) 2 ||Co fc Jo(ji fc r)||| 2( 



(6.6) 



- II ,°ll 2 
fc=i 



and convergence in the vanishing viscosity limit follows from the condition 
in Equation (|2.2p of Theorem 12.31 (Examining the argument in |6] shows 
that the convergence rate in Equation (|2.ip is bounded by C(z^t) 1 / 2 .) □ 

We could have concluded the proof another, more indirect way, as follows. 
From Equation (16. 5p we see that uj(t) is continuous in the L 2 (0)-norm at 
time zero. Thus, we can multiply Equation (|6.4f) by uj and integrate over 
time to give 



(n) + v ||Vw(t)||^ (n) = "J r (^ ■ n)uj. (6.7) 
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But = —(An) = on r because An is in V as we showed, and 
we conclude by integrating over time that u> is in L°°([0, oo); L 2 (Q)), and 
Equation (12. ip follows from the condition in Equation (|2.2p of Theorem 12.31 

This argument shows that whenever the vorticity is continuous at time 
zero in the L 2 (f2)-norm and Au(t) lies in V for all t > 0, Equation (|2.1j) 
holds. The latter condition, however, is very special. In fact, for n° with 
regularity as in Theorem 16.21 dtu will always be in V, so An = (dtu + n • 
Vn + Vp)/v is in V if and only if n • Vn + Vp is in V. This, in turn, can 
hold only if Vp = — n • Vn on T. 

This suggests that for initial velocities for which u(t) is continuous in 
the L 2 (f2)-norm at time zero we might attempt to make an argument using 
Equation (|6.7p . though now without the right-hand side vanishing. This 
would require either control on Vw, as we had above, or on u itself. (For 
solutions to (NS) with the boundary condition oj = u ■ n = rather than 
n = this term would vanish and one can obtain the vanishing viscosity 
limit easily — though not using Theorem 12.31 which one would need to show 
applies to such solutions — as done, for instance, in and |12j.) 

Another line of attack is also suggested by the proof of Theorem 16.21 — 
to gain control on the coefficients g(t) either by assumptions on the initial 
velocity or on the solution. Unless bounds as remarkably strong as those 
obtained in the proof of Theorem 16.21 are achieved, though, something more 
sophisticated must be employed to obtain Equation (|2.ip . 

7. Interpretation in terms of length scales 

In P|, Cheng and Wang consider the vanishing viscosity limit in the setting 
of a two-dimensional rectangular channel R, periodic in the x direction with 
period L and with characteristic boundary conditions (which include no-slip 
boundary conditions as a special case). They decompose any vector n on 
R of sufficient regularity as n = ^JLo e 2m '^ x / L u^ and define the projection 

p ku = Ej=o onto the space spanned by the first k modes. This 

in effect allows one to isolate successively finer-scale spatial variations in the 
direction tangential to the boundary. They then construct an approxima- 
tion sequence {v L } to n by letting v L be the solution to the equation that 
results after projecting each term in (NS) using P]y. (We have changed 
their notation somewhat.) Their v is the approximate-solution analog of 
the exact solution truncation represented by u L in Equation (|2.9p . 

The main result in [3] is that v L ^ v ' converges to u in L°°([0, T]; L 2 (S7)) 
as v — * 0. The requirement on L(y) is the same as our condition on L in 
Equation (|2.10p (with the additional condition that L(v) — > oo as v —> as 
one would expect), so convergence of v L ^ to u occurs when only tangen- 
tial length scales of order larger than v are included in the approximations. 
(All length scales in the normal direction, however, are included. See Re- 
mark (|5.ip concerning this issue in regards to the vorticity.) 
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The result in [3] makes an important observation about the difficulty of 
determining numerically whether or not the vanishing viscosity limit holds. 
Our method of decomposing the solution using the eigenfunctions of the 
Stokes operator, on the other hand, says little about computation, since 
approximating this decomposition numerically is probably as least as hard as 
approximating the solution itself. Nonetheless, it more directly characterizes 
the properties of the solution itself at different length scales. 

The analog to the result in [3] is Theorem 15.51 which shows that Temam 
and Wang's conditions in Equation (12. 5ft and Equation (|2.6p . when applied 
only to the modes with tangential wavelengths of CL(v) or higher, holds 
as long as the condition on N in Equation (|2,10p hold. This does not, 
however, imply that it^M converges to u in the vanishing viscosity limit, 
only that if the vanishing viscosity limit fails to hold, the failure originates 
in the behavior of the tangential component of the gradient projected into 
the space spanned by the modes with tangential frequencies of order L(y) 
or higher; that is, at length scales of order v or lower. 

The other conditions in Theorem 12.41 give alternative ways to measure 
the behavior of the solution at different length scales or frequencies. They 
show that we cannot simply say that if the vanishing viscosity fails to hold 
then the failure lies in the behavior of the solution at any particular range 
of length scales, but rather that the pertinent range of length scales varies 
with the measure of behavior. Whether any of these conditions brings us 
any closer to proving that the vanishing viscosity limit holds in general for 
smooth initial data in a bounded domain or to proving that it fails to hold 
in at least one instance remains completely unclear. 

Appendix A. Bounds on the Eigenfunctions 

In Lemma IA.1I we state the basic identities involving the Bessel functions 
that we use. We then give a series of lemmas that lead to the bounds on 
the velocity and vorticity of the eigenfunctions in the boundary layer that 
we used in the proof of Theorem 12.41 

It is perhaps important to note that in the proofs that follow we avoid 
the use of asymptotic formulas for the Bessel functions, even when such 
formulas might appear to be useful. This is because we need to deal with 
the relative values of Bessel functions of different orders near a zero of one 
of the Bessel functions, and it is precisely in these situations that the errors 
in the asymptotic formulas dominate. Also, most of the following lemmas 
apply without change to their proofs with n being any nonnegative real 
value. 

Lemma A.l. For all nonnegative real numbers n and x, 

2nJ n (x) - xJ n _i(x) = xJ n+ i(x), (A.l) 



2J' n (x) = J n -i(x) - J n+1 (x), 



(A.2) 
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Jn-l{x) = ~Jn{x) + J' n {x), (A.3) 

J n+1 (x) = -J n {x) - J' n {x), (A.4) 
x 



x n J n (ax) 



a 



J x n J n -i(ax) dx, (A. 5) 



xJ n (ax) 2 dx = — 



x 2 J' n (ax) 2 + ( x 2 j ) Jn( ax Y 



(Ai 



J n (ax)x n = —a J J n+ i(ax)x n dx, (A. 6) 

(J3 2 — a 2 ) I xJ n {ax)J n {[3x) dx = x \aJ' n {ax)J n {[3x) — (3J' n (f3x)J n (ax)^ , 

(A.7) 

a 2 

f x 2 

/ xJ n (ax) 2 dx = — [J n (ax) 2 - J n _ 1 (ax)J n+ i(ax)] . (A.9) 

Proof. These are standard identities for Bessel functions. For instance, see 
Equations (6.28), (6.29), (6.30), (6.31), (6.38), (6.39), (6.51), (6.52), and 
(6.53) of [2]. □ 

Lemma A. 2. For all nonnegative integers n and all positive integers k, 

1 < jn+l,k - jnk < ^! 

where j n k is defined in Equation |^.^[). 



Proof. Let j u k be the k-th positive zero of J v , where we now allow v to be 
a real number in the interval [0, oo). It is shown in [S] and [2] that for all 
k > 1, j u k is strictly concave as a function of v and that dj u k/dv > 1 (see 
also [5]). Thus, the function n i— > j n+ \^ — jnk is strictly decreasing as a 
function of n. But by Equation (2.9) of [3], j n +\ : k — jnk — > 1 as n — > oo, so 

jn+l,k — jnk > 1- 

The positive zeros of Jo lie in the intervals {mir + |-7r,m7r + |vr), m = 
0, 1, . . . , and the positive zeros of J\ lie in the intervals (ra'ir + ^ir, m'ir + 
^7r), m' = 1,2, ... . That the zeros lie in only these intervals is shown in 
Section 15.32 p. 489 and Section 15.34 p. 491 of [18] using an approach 
of Schafheitlin's. That each of these intervals contains at least one zero is 
shown on p. 104 of [2] . But j\k —jok > 1 as we showed above so each interval 
contains precisely one zero. Because the zeros of Jo and J± are interleaved 
(see p. 106 of [2], for instance) we can then conclude that jik — jok < § ■ But 
as we observed above, the function n i— > j n +i,k — jnk is strictly decreasing as 
a function of n, so j n +i,k — jnk < f holds for all n > 0. □ 
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Lemma A. 3. For all n = 0, 1, 



and k = 1, 2, . . . , 



n + < j n k < 7r(n/2 + fc) < 7r(n + A;). 



Proof. By Lemma lA.2( for all n and j, 

n 

infc = Jok + X] C?mA - 3m-l,k) > Jok + n> n + k, 



m=l 



because jok > k (which follows directly from Equation (14, 3p : see p. 485-486 
of [18], for instance). By an observation in the proof of Lemma [A.2l it follows 
that jok < irk, and a similar argument using the inequality j n +i t k — jnk < § 
from Lemma I A. 21 gives the upper bound on j nk . □ 



Lemma A. 4. Let a 

1,2,..., 



j n +i,k and (3 = j nk . For n = 0, 1, 2, . . . and k 



J n {ax) 



J n (a) 





< 1 if - < x < 1. 
a 



Proof. Let g(x) = J n {ax)/\ J n (a)\. From Equation (jA.4j) . J' n {a) = (n/a)J n (ot) 
so J' n {oi) has the same sign as J n {a). From this we conclude that \g\ is in- 
creasing in a left-neighborhood 7V of 1. 

Between each zero of J n there is exactly one zero of J n +i (see p. 106 of 
[2], for instance). Between each zero of J n there is also exactly one zero of 
J' n , because the maximum values of J n are all positive and the minimum 
values are all negative (see, for instance, p. 107 of [2]) and J' n has no repeated 
positive roots (this follows from the defining equation Equation (|4.3p ). Thus, 
the neighborhood N includes all x such that j3 < ax < a. Since \g(l)\ = 1 
it follows that \g(x)\ < 1 for all such x. □ 



Lemma A. 5. Let a 

that for all n = 0,1, . 

J n+ i(ax) 



Jn+i.fe and (3 = j nk . 
and k = 1, 2 . . . n, 



Jn(a) 



< Cn(l - x) if — < x < 1 



There exists a constant C such 



a 



Proof. Since J n+ \{a) = 0, Equation (|A.5[) with n + 1 in place of n gives 



J n+ i(ax) 



a 



r n+l 



t n+i J n {at)dt. 



As long as f3 < ax < a, J n (at) does not change sign on the interval (x, 1] 
and has its maximum value on this interval at 1, as observed in the proof of 
Lemma IA.4I Thus, 



|J n+1 (ax)| < 



a 



r.n+1 



\J n (a) 



4.n+l 



dt < 



a\J n (a 



\n+l v ' 



But by Lemma IA.2I and Lemma IA.31 

/3 < vr/2 



a 



a 



a 



n + 2 



P 



a 



> 1 



7T 



2n + 4 
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SO 



(/3/a)~ (n+1) < f 1 



7T 



>+l) 



< e^ 2 . 



2n + A / 

the last inequality following from elementary calculus. We conclude that 

\J n +i(ax)\ < Cn \J n (a)\ (1 - x), (A.10) 

which completes the proof. □ 

Lemma A. 6. Let a = j n +ik a nd (3 = j n k- There exists a constant C such 
that for all n = 0, 1, . . . and k = 1,2 ... ,n, 



J n -i(ax) 



J n (a) 



< C if - < x < 1. 
a 



Proof. Because the positive zeros of J n -\ are interlaced with those of J n , 
J n -i does not change sign on the interval [/3, a]. From Equation (IA.4j) with 
n — 1 in place of n, J/ l _ 1 (/3) = ((n — I)/ (3)J n _i((3), so J' n _i{(3) has the same 
sign as J n _i(/3), and we conclude that J n -\ reaches it maximum value on 
the interval \p,a] at j3. Therefore, for (3/ a < x < 1, 



J n -i(ax) 



J n (a) 



< 



Jn-l(P) 



J n {a) 



But, by Equation {AT]), J n +l((3) = 2(n//3) J n (P) - J n -l{0) = -Jn-l(P), so 



J n ~i(ax) 



J n (a) 



< 



J n (a) 



< Cn(l - (3/a) < C, 



where we used Lemma IA.5I and Lemma IA.3I 



□ 



Lemma A. 7. We have ||o>i 



< 25 when 5 < A 

1/2 



-1/2 



Proof. Let tuj = u) nk and a = j n +i,k = X/ ■ Then 



1-6 



rJ n (ar) 2 dr 



1-5 



' Jn(a) 2 



dr. 



^Il£2(r 5 ) 

In the integrals above, with (3 = j n k, 

P/a = l-(fi-a)/a<l-l/a = l- XJ 1/2 <l-5 <r <1, 

where we used Lemma lA.2} and the lemma follows from Lemma IA.4I □ 

Employing Lemma lA.71 we can extend its range of applicability, though 
with a higher bound on the width of the boundary layer. 

— 1/2 

Lemma A. 8. For all n = 0,1, ... , k = 1, . . . , n, and all 5 < 2irX nl , 

1 1 2 

ll w nfc|lL 2 (r 5 ) — 

Proof. It follows from Lemma lA.31 that j n ,n/jn,i < + n)/(n + 1) < 2ir; 
the lemma follows from this inequality and Lemma lA.71 □ 
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Remark A.l. It is possible to extend Lemma IA.8I to include all values of 
k. The idea of the proof is that for k > n, u> n k passes through k complete 
half-periods (annuli in the unit disk lying between successive nonnegative 
zeroes of JnUn+i,^)) and ends with a partial period. Since J n (x) decays like 
x 1 / 2 and the spacing between consecutive zeros of J n approaches a constant, 
the L 2 -norms of u> n k on each of those half-periods converges to a constant, 
and since the L 2 -norm of u> n k on the entire unit disk is 1, the square of the 
L 2 -norm of Lu n k on the last half-period is less than C/k (with C near 1). 
But the last half-period has a width that is greater than C/k. Extending 
this argument to m periods, what we have shown is that 

II^HWcw*) " m/k - 
With the assumed bound on 5, we choose m so that m/k is of the same 
order as 5, and the proof is essentially complete. 

Lemma A. 9. There exist positive constants C\ and C% with C2 < 1 such 
that for all n = 0, 1, . . . and all k = 1, . . . , n, 

1 1 2 3 
Ipn&llizg^) < C\5 

when 5 < C^A,^ 2 . 

Proof. In the proof that follows, we will often use Lemma lA.31 without ex- 
plicit mention. 

Let a = j n +i,k- We bound first the radial component of u n k- We have, 
J n {ar) - J n {a)r n 



J n (a)r 



_1 5n(r), 



where 

= J -l 

J n (a) J n (a) 



J n (ar)r n a f 1 J n+ i(ax) f 1 B nk (x) 
9n{r) = —— 1 = — — - / dx = -a — dx > 



and 



B n k{x) 



J n+ i(ax) 



J n (a) 

To verify the second equality we use the identity in Equation (IA.6jl . from 
which it follows that the second expression for g n is 

J n (ar)r~ n + C 
J n {a) 

for some constant C. But all three expressions for g n are zero at r = 1, so 
we have the correct limits of integration in the second expression. It follows 
from Lemma IA.5I and our third expression for g n that 

f 1 (In 2 
\9n(r)\<CnaJ x" n (l - x) dx < ^-(1 - r) 2 

for all 1 — r < 1/a. 
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From Equation (14,81) . 



g(r)r 



n—l 



so when 5 < 1/a we have 

KfcllL2(r 4 ) = 27r / r Kkl ^ < — - / „2n-l * 

< Cn 2 {l - d) 1 " 2 ™^ < C5 3 . 
In the last inequality we used 

<5 < 6 2 A nl = < 



3n+l,k n + 1 
so 



(l_jfn-i> (1-^-1 = (G(C 2 ,n + l)r(l 



^ \ 2n— 1 / /-y 

Co \ ,„,„ . sv9. / . Co 



n+1 J \ n + 1 

> (1 - C 2 ) 2 (1 - C 2 )" 3 = (1 - C,)" 1 = C7 > 0, 

where G is the function of Lemma IA.11I 

For the angular component of u n k, we write 

aJ n+ i(ar) - aJ n _i(ar) + 2nJ n (a)r n ~ 1 

= [2nJ n (ar) — arJ n ^i(ar)] + arJ n ^\{ar) — 2nJ n {ar) 

+ aJ n+ i(ar) - aJ n -i(ar) + 2nJ n (a)r n ~ 1 

= a(r + l)J n+ i(ar) + 2n [^(a)^" 1 - J n {ar)} 

+ aJ„_i(ar)(r - 1). 

From Equation (14. 8p we then have 



a 4 Jn{a) 2 a 4 \ J n {oi) 

+ C 



a 2 Jn-i(ar) 2 (l-r) 2 ( A11 ) 



a 4 J n (a) 2 

<C(l-r) 2 + -^(r"-V-i(r)) 2 , 

where we applied both Lemma I A. 51 and Lemma lA. 61 
The first term in Equation (jA.lip contributes no more 



C r(l - r) z dr<C (1 - r) 2 dr < C8 6 , 
Jl-S Jl-5 

and the same is true of the second term in Equation (jA.lip arguing as for 

u r nk , and this completes the proof. □ 

Lemma A. 10. When m^n, {u mj ,u nk ) = (uJmj,^nk) L 2(r c ^ = 0. 
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Proof. We have, 

{u: mj ^ nk ) L2{Tcv) = J 1 r/(r) d6 dr, 

where f(r) is a product of two Bessel functions. When m ^ n, the inner 
integral is zero. A similar argument gives (u m j, Unk) jjz^y 

Lemma A. 11. Xei a be in (0, 1) and define G a : [1, oo) — > [0, oo) fry 

C.(,)-(l-f)". 

T/ien /or aZZ x > 1, 

1 - a < G a (x) < e _a . 
Proof. The proof is elementary. □ 
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